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We discuss topological objects, in particular the non-Abrikosov vortex and the magnetic knot 
made of the twisted non-Abrikosov vortex, in two-gap superconductor. We show that there are 
two types of non-Abrikosov vortex in Ginzburg-Landau theory of two-gap superconductor, the D- 
type which has no concentration of the condensate at the core and the N-type which has a non- 
trivial profile of the condensate at the core, under a wide class of realistic interaction potential. 
Furthermore, we show that we can construct a stable magnetic knot by twisting the non-Abrikosov 
vortex and connecting two periodic ends together, whose knot topology iT3(S 2 ) is described by the 
Chern-Simon index of the electromagnetic potential. We discuss how these topological objects can 
be constructed in MgB2 or in liquid metallic hydrogen. 
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I. INTRODUCTION 



Topological objects, in particular finite energy topo- 
logical objects (monopoles, vortices, skyrmions, and 
knots), have played increasingly important role in physics 
[H 13- LJ- CO- 0] ■ In condensed matter the best known topo- 
logical objects are the Abrikosov vortex in one-gap super- 
conductors and similar ones in Bose-Einstein condensates 
and superfluids, which have been the subject of intensive 
studies. A recent advent of two-component Bose-Einstein 
condensates and two-gap superconductors 0, Q, how- 
ever, has opened up an exciting new possibility for us to 
construct far more interesting topological objects in lab- 
oratories. It has already been shown that non-Abrikosov 
vortices whose topology is fixed by ^{S 2 ) and finite en- 
ergy topological knots whose topology is fixed by tt z(S 2 ) 
exist in these condensed matters USEES HQ- The 
reason for this is that these condensed matters are made 
of two components which can be viewed as an SU (2) 
multiplet. In general this type of topological objects is 
possible when one has a multi-component condensates, 
which allows the non-Abelian topology. 

The purpose of this paper is to discuss new topological 
objects in Ginzburg-Landau theory of two-gap supercon- 
ductor in detail. With a most general U(l) x U(l) sym- 
metric potential which can describe a wide class of two- 
gap superconductors we first show that there are two types 
of non-Abrikosov vortex, D-type and N-type, in two-gap 
superconductor. The D-type has no concentration of the 
condensate at the core, but the N-type has a non-trivial 
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profile of the condensate at the core. The reason why the 
two-gap superconductor has two types of vortex is that 
the vortex in two-gap superconductor allows two differ- 
ent boundary conditions. In terms of topology the non- 
Abrikosov vortex is described by two types of topology, 
non-Abelian ^{S 2 ) topology or Abelian ^(S 1 ) topology. 
And within the same topology both D-type and N-type vor- 
tices exist. In particular, there are infinitely many D-type 
vortices classified by the natural number k which have the 
same topology. Moreover, the magnetic flux of these non- 
Abrikosov vortices can be integral or fractional, and the 
integral flux vortex has the tt2(S 2 ) topology and the frac- 
tional flux vortex has the tti(S 1 ) topology. We show that 
the N-type vortex has a 27rn/<7-flux or a fractional flux 
(a fraction of 2irn/g), but the D-type vortex has 2irk/g 
more flux than the N-type vortex. These characteristic 
features of the non-Abrikosov vortex are clearly absent 
in the Abrikosov vortex which carries 2nn/g-&xx. whose 
topology is fixed by 7r 1 (S' 1 ). 

Next, we show that the non-Abrikosov vortex can be 
twisted to form a helical vortex which is periodic in z- 
axis. More importantly, we show that we can construct 
a stable magnetic knot in two-gap superconductors by 
smoothly bending the helical vortex and connecting the 
periodic ends together. The vortex ring acquires the knot 
topology iT3(S 2 ) which is fixed by the Chern-Simon index 
of the electromagnetic potential. Because of the helical 
structure of the magnetic flux the knot has two magnetic 
flux linked together, one around the knot tube and one 
along the knot, whose linking number is given by the 
knot quantum number. And the flux trapped inside the 
vortex ring provides a stabilizing repulsive force which 
prevents the collapse of the knot, because it can not be 
squeezed out. This means that the knot has dynamical 
(as well as topological) stability. 
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It is well-known that multi-gap superconductor may 
have interband Josephson interaction [14| • We consider a 
most general quartic Josephson interaction in two-gap su- 
perconductor, and show that the presence of the Joseph- 
son interaction does not affect the existence of the above 
topological objects, but can alter the shape of the solu- 
tions drastically. We show that in the presence of the 
Josephson interaction we have a magnetic vortex which 
can be viewed as a bound state of two fluxes, which be- 
comes a braided magnetic vortex when twisted. 

The paper is organized as follows. In Section II we dis- 
cuss a most general quartic potential in Ginzburg-Landau 
theory of two-gap superconductor in mean field approxi- 
mation, and study the vacuum structure. In Section III 
we show that the Ginzburg-Landau theory of two-gap su- 
perconductor can be understood as a theory of CP 1 field 
coupled to a scalar field and the electromagnetic field, 
and argue that the topology of the theory can be de- 
scribed by the CP 1 field and the electromagnetic field. In 
Section IV we construct the non-Abrikosov vortex in two- 
gap superconductor, and show that there are two types 
of boundary condition which allow two types of magnetic 
vortex, the D-type which has no concentration of conden- 
sate at the core and the N-type which has a non-trivial 
concentration of condensate at the core. Moreover we 
show that there are two types of topology, non-Abelian 
TT2(S 2 ) and Abelian 7Ti(5 1 ), which describes these vor- 
tices. We show that the magnetic flux of these vortices 
can be integral or fractional depending on the parameters 
of the potential, but the D-type vortex has 2i:k/g more 
flux than the N-type vortex. In Section V we show that 
we can construct a helical magnetic vortex in two-gap su- 
perconductor, by twisting the non-Abrikosov vortex and 
making it periodic in z-axis. In Section VI we construct 
the magnetic knot bending the helical vortex and con- 
necting the periodic ends together, and show that the 
knot topology n^(S 2 ) is described by the Chern-Simon 
index of the electromagnetic potential. In Section VII 
we consider the Josephson interaction, and show that the 
inclusion of the Josephson interaction does not affect the 
existence of the topological objects in two-gap supercon- 
ductor but alter the shape of the solutions drastically. In 
Section VIII we discuss the non-Abelian superconductiv- 
ity which can describe a two-gap superconductor made of 
two condensates which carry opposite charge, and argue 
that the non-Abelian superconductivity can be realized 
in liquid metallic hydrogen (LMH). Finally in Section IX 
we discuss the physical implications of our results, and 
discuss how one can identify these topological objects in 
MgB 2 and LMH. 



II. EFFECTIVE POTENTIAL OF TWO-GAP 
SUPERCONDUCTOR 



In mean field approximation the free energy of the 
two-gap superconductor could be expressed by 0, Q, 
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where V is the effective potential. We choose the poten- 
tial to be the most general quartic potential which has 
the U(l) x U(l) symmetry, 
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where Ay are the quartic coupling constants and fii are 
the chemical potentials of 0j (i — 1,2). One might like to 
include the Josephson interaction to the potential which 
breaks the U(l) x U{1) symmetry down to U(l). The 
Josephson interaction will be discussed separately in the 
following. But as we will see, the inclusion of the Joseph- 
son interaction does not alter the qualitative features of 
the topological objects we discuss in this paper. 
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where V is the normalized potential 
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From now on we will assume that all coupling constants 
except A12 are positive. 

To find the vacuum of the potential let 



dV 
dV 



An|0i| 2 + A 12 |^ 2 | 2 — Mi = 0, 
M2 = 0, 



Ai 2 |0i| 2 + A 22 |02| 2 



and find the extremum 

|0i| 2 =0 2 , 



I0 2 | 2 
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To check whether this extremum is the maximum or min- 
imum, consider the Hessian 

d 2 V 2 
dCtg = det ^P^|^ AllA22 " A ' 2 

= A. (9) 

There are three possibilities; positive, zero, or negative 
A. We consider each case separately. 

A. A > 0: In this case we have \\ 2 < A11A22, and the 
extremum © becomes the local minimum. But since \<f>i\ 
have to be positive we have the following vacuum 



< |0i| > 

< \<h\ > 




for A12 < or for 



0<A 12 , £l<£<£i. 

A22 M2 A12 
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Notice that both < \<j>i\ > and < |</> 2 | > are non- 
vanishing. But for 



0<A 12 , ^<^<tl. 

A22 A12 M2 

we have the following vacuum from © , 



< I0i| >) = 

< \<h\ > ) I 



(12) 
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Finally, when 



U < Ax2, — < t — < t — , 
M2 A22 A12 
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we can always transform this case to the case (1 1 121) by re- 
labeling 4>i and (f>2 as $2 and <f>i, so that in this case we 
can assume that the vacuum is still given by (|13fl without 
loss of generality. 

B. A = 0: In this case we have A 2 2 = A11A22, and the 
potential JSJ) is reduced to 
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So for A12 < 0, the potential becomes unbounded from 
below, so that it has no minimum. For 



U < A12, - — — - — < — , 

A12 A 2 2 fl2 



we have the following vacuum 



< I0il > ) = / VMi/AiT 

< \<h\ > } \ 



Next, consider the case 



0<A 12 , ^<^i = ^i. 

(12 A12 A22 



(16) 
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But this can be transformed to (|16fl with the re-labeling. 
So we can assume that the vacuum is given by l|17f) with- 
out loss of generality. Finally, when 



Ml _ An _ A 12 
M2 A12 A22 

we have the degenerate vacuum 

Mi < l^i I > 2 +M2 < 1 02 1 > 



2_ MlM2 
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This case includes the special (and familiar) SU(2) sym- 
metric case 



An = A12 = A 22 = A, u\ = M2 = M ; 
<|0i|> 2 + <|0 2 |> 2 =^. 
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In this case the Hamiltonian Q has the full SU (2) sym- 
metry. 

C. A < 0: In this case we have A 2 2 > A11A22, and 
the extremum (|SJ) becomes the local maximum. So the 
minimum state must satisfy 



l0l| 2 |0 2 | 2 =O. 

(14) Now, by inspection one can show that when 



0<Ai2, ^< 
A12 



' An Mi 
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the vacuum must be 



< I0i| > 1 = ( VMi/An 

< I0 2 | > 



Next, consider the case 
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This case can be reduced to the above case (by re-labeling 
(15) <pi and <p 2 ), so that when A12 is positive one can assume 
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that the vacuum is given by (|24|l without loss of gen- 
erality. Finally when A12 is negative the potential has 
no minimum, because it is unbounded from below. This 
must be clear from (I15|l . 

In summary, we have three types of vacuum state: 
A. Type I: Integer flux vacuum 



< I0i| > 



(26) 



This is possible when we have one of the following three 
cases, 



(a) < A12, 

0) < A12, 
(c) < Ai 2 , 



A 12 -Vu < Mi 

A22 A12 /i2 
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We call this integer flux vacuum because, as we will see, 
for this type of vacuum the magnetic vortex has an inte- 
ger flux. 

B. Type II: Fractional flux vacuum 
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<\h\> J V02 
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This is possible when we have one of the following three 
cases, 



(a) A12 < 0, 
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(c) A12 = 0. 
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We call this fractional flux vacuum because, as we will 
see, for this type of vacuum the magnetic vortex has a 
fractional flux. 

C. Type III: Degenerate vacuum 



Hi < |0i I > 2 +M2 < 1 02 1 > 2 = 



M1M2 
A12 



This is what we have when 
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Notice that the potential © has no vacuum when 
Al2 <0, iiL<l Al2 l 



IA12I 



A2 
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All other cases can be reduced to one of the above cases 
by re-labelling 0i and 02. As we will see the vacuum 
structure will play an important role in the following. 



Notice that with 



An + A22 + 2Ai 2 
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(33) 



we have 



An = A + (3 + a, A 2 2 = A + /3-a, 
Ai2 = A-/3, 
Mi=M + 7, M2=M~7, (34) 

so that the potential can be written as 
2^' ~ A ' 2' 



^(l0l| 2 + l02| 2 -^) 2 + ?(|0l| 4 -|02| 4 ) 
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In terms of the new parameters wc have 



2(Au + 7 A)-a(M + 7) 



2(/3m - 7A) + a(M - 7) 
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where A = 4/3 A — a 2 



III. DYNAMICS OF TWO-GAP CONDENSATES 



With this preliminary one may study the topological 
objects of two-gap superconductor minimizing the free 
energy. On the other hand, to study a static solution, 
one might as well start from the following relativistic 
Ginzburg-Landau Lagrangian 



£ = -\D 



= (dp + igAfJcf), 



(37) 



which reproduces the the free energy J3J in the static 
limit. The Lagrangian has the equation of motion 



D^i = 



D 2 (j>2 = 
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To understand the meaning of this we let 

0=-j|P£, I0 f 0l = f, ^=1, 
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and find the following identities 



we can express the Hamiltonian @ as 



(9,n) 2 =4(|^| 2 -|^| 2 ), 

1 2i 
—h ■ (d^n x d u h) = —(d^d u £ - d^d^) 



9 

= d^Cn — d v C^, 
\{igCn- a ■ d^h) 



c„ 



(40) 



From these we can reduce (|3*8"|) to 

2 P- (\{d,nf + g\A,- l -C,) 2 )p 



h x <9 2 7i + 2^^-h x 9 u n 



d il F^=j v = g*p i (A v -~C 1 



(41) 



This is the equation for two-gap superconductor, which 
allows a large class of interesting topological objects, 
straight magnetic vortex, helical magnetic vortex, and 
magnetic knot, all with Air / g-fhxx, 2w / g-ftux., or fractional 
flux. 

The equation i|38|) is an equation of the complex dou- 
blet 4> which has four degrees. But notice that the equa- 
tion (|41|l is, except for C^, expressed completely in terms 
of the CP 1 field n and the scalar field p. Moreover, 
(I40JI tells that C M can also be written in terms of h. In 
fact n uniquely defines a righthanded orthonormal frame 
(ni,fi2,n), with h% X rii = h, up to the U(l) rotation 
which leaves h invariant. Then is given (up to a 
U(l) gauge transformation) by the Mermin-Ho relation 



i 



— hi ■ 8^2, 
9 

d fl C u -d u C IJ , = --h-(d fJl hxd u h). (42) 

This tells that we can transform the equation l|38fl of the 
complex doublet condensate <fi to the equation l|41|l of the 
CP 1 field h and the scalar field p. In fact, with 



n = -(d^pf + -g 2 p 2 Bl + -p 2 (9» 2 + V(p,n 3 ) 



where 



1 1 

+ ^ [G/xk - — h ■ (d^h x d v h)] ' 



V(p,n 3 )^^(l + ^n 3 + ^nl)p^ 
--(l + -n 3 )p, 



(44) 



(45) 



and n 3 — £^£1 — This means that the Ginzburg- 

Landau theory of two-gap superconductor can be under- 
stood as a theory of CP 1 field h (coupled to p and B^) 
This is because the U{1) gauge invariance 
of (f3*7| reduces the physical degrees of the complex dou- 
blet (j) to p and h, and the massive photon B^. As we 
will see, this has a very important physical implication, 
because this tells that the topology of two-gap supercon- 
ductor can be described by the topology of h and A M . 

The equation (|4 1|> allows two conserved currents, the 
electromagnetic current and the neutral current fc M 

El, 



k, = gipi [A^Ui - 86) + -fa&ti 



(46) 



which are nothing but the Noether currents of the U(l) x 
U(l) symmetry of the Hamiltonian Q. Indeed they are 
the sum and difference of two electromagnetic currents 
of <j>i and <f>2 



k - 



J (2) - 



(47) 



Clearly the conservation of follows from the last equa- 
tion of (|4ip. But the conservation of k^ comes from 
the second equation of (|41|l . which (together with the 
last equation) tells the existence of a partially conserved 
SU(2) current % O, 



„ /l 

.9P y-^n x 9 M n 



(48) 



This SU (2) current is exactly conserved when a — (3 — 
7 = 0. But notice that fc M = k ■ j^. This assures that 
we have the conservation of k^ even when a(3"f 0. It 
is interesting to notice that and fc M are precisely the h 
and k components of j L . 



5,, 



G>„ = d^Bu — d v B^, 
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IV. NON-ABRIKOSOV VORTEX 



we can express (|41|l completely in terms of n, p, and _/3 M . 
This is not accidental. Indeed with JSSJl, gDJl, and |@3J|, 



The two-gap superconductor allows different types of 
interesting magnetic vortex || 0, 0, . In terms of the 
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structure of the vortex they are classified to two types, 
the D-type vortex which has no concentration of the con- 
densates at the core and the N-type which has a non- 
vanishing concentration of the condensates at the core. 
The reason for this is that, unlike the Abrikosov vortex 
in ordinary superconductor, the vortex in two-gap super- 
conductor allows two different boundary conditions at 
the core. 

To discuss the straight vortex let (g, ip, z) be the cylin- 
drical coordinates and choose the ansatz 



pie), £ = 



/ fie) , ■ x 
cos — - — exp(— vnipj 



V 



2 . fie) 

sin — - — 



A„ = -A(g)d^. 



With this we have 



(sin /(g) cos nip 
sin /(g) sin nip 
cos /(p) 

cos /(g) + 1 

Cf, = n d^ip, 
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2 ( , COS / + 1 \ 

= n 9P A df.ip, 



k v = ngp 2 



{A cos / 



cos / + 1 



du<P, 



and the Hamiltonian (|44|l becomes 



n2 P 2 ( A cos/ + l N| 2 + n 2 ^, 



2g 2 g 2 
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(P - P ) + jiP - —)P cos/ 



4p 4 cos 2 / 



/£ 

2A 



(51) 



With this lHJl becomes 
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(p 2 -p 2 ) + ?(p 2 -^)cos/ + fp 2 cos 2 / 
A a A 



(>■ 



A 



(2 7 -(| +/ 3cos/V 2 ) sin/, 

1 : 2 2 / a COS / + 1 

A-ffV(4 1 



>J 



0. 



(52) 



Notice that this can also be derived by minimizing the 
Hamiltonian l|51|) . 

To solve the equation, we have to fix the boundary 
conditions. To determine the possible boundary condi- 
tion at the core we expand p(g), /(g), and A(g) near the 
origin as 



p(g) ~ po + pig + p 2 g 2 + pzg 3 + 
/(g)~/ + /ig + / 2 g 2 + /3g 3 + 
A(g) ~ a + dig + a 2 g 2 + a 3 g 3 + 



and find that the smoothness at the core requires 



(53) 



^((2a - cos/o - l) 2 + sin 2 / ) + p (^(2a - l)/ a sin/ + (2a - cos/ - l)a x ) 

h^((2a - cos f - l) 2 + sin 2 / - 4)1 + f- + % ((2a„ - cos f - l) 2 + sin 2 f - g fc - 2 
4V n^/J g L 4 V n z /J 



^((2 ao -cos/ -l) 2 + sin 2 / 

(2an — l)sinfn „ 1 N \ 1 

V g2 ; + (((2ao - 1) cos/o - ^)/i + 2oi sin/oj - + ... = 0, p ^ 

(2a n — l)sinf n 3, \1 

^-5 ^ ^ + 2a - 1 cos /o - — /i + 2ai sin/ - + ... = 0, Po = 0, Pl ^0 

g l V v n z/ / g 



0. 



(2a - l)sin/ 



((2oo - 1) cos/o - 2fe "a" 1 )/i + 2aisin/ )- 



= 0, po = Pi = —Pk-i = 0, pk ^ 



oi , 9 2 Po 



+ ( 2a o - cos / - l) - (3a 3 - g 2 p p 1 (2a - cos/ - 1) 
-5 2 Po(«i + ^/isin/ ))g+ ... = 0. 



(54) 



7 



Now, consider the case po ^ first. In this case the 
last equation requires a\ = and 2ao — 1 = cos/o, and 
the first equation requires sin/o = and p\ = 0. With 
sin /o = we must either have /o = tt and a = or 
/o = and oo = 1- And with /o = 7r, we have /i = 
for n ^ ±1 from the second equation. One might choose 
fo = instead of / = tt, but this does not lead us to 
a new solution. Next, consider the case po = 0, pi 7^ 0. 
In this case the first two sets of equations tells that we 
may have fo — n, ao = ±l/n, and /1 = for n 7^ =pl or 
n 7^ ±3. Similarly, for p = pi = ... = p k -i =0, pk ^ 0, 
H54JI tells that we may have fo = 7r, ao = ±fc/n, and 
fx = for n ^ Tl or n ^ ±(2fe + 1). This tells that we 
can choose the following boundary condition at the core 
for the vortex described by the ansatz (|49|l in two-gap 
superconductor [1,0,0: 
A. Dirichlet boundary condition 

p(0) = 0, p(0)^0, 
A(0) = ~, /(0) = tt, f = forn^l. (55) 

In particular, for n = 1 we must have A(0) = — 1. In 
general, the Dirichlet boundary condition can be written 



t(o) = o, 



do' 



fc (o)^o, 



A(0) = - 



k 



dg k 

/(0) = tt, / = 0/orn^l, (56) 



so that for n = 1 we must have A(0) = — fc. 
B. Neumann boundary condition 

p(0)^0, p(0)=0, 
A(0)=0, /(0)=7r, f = 0forn^l. (57) 

So here we have A(0) = for all n. This shows that 
the magnetic vortex in two-gap superconductor allows 
two types of boundary condition which arc different from 
what we have in ordinary superconductor. This is a new 
feature of two-gap superconductor which will have a deep 
impact in the following. 

To determine the boundary condition at the infinity 
notice that at the infinity all fields must assume the vac- 
uum values. In particular, the electromagnetic current 
must vanish at the infinity. This means that we must 
have 

p(oo) = V2(< [^1 > 2 + < \<h\ > 2 ), 



cos/(oo) =< n 3 >= 



< 01 >' 



< 01 >< 



< 02 >< 



A(oo) 



«/(c 



<l0l|> 2 



<|0l|> 2 + <|0 2 |> 2 



(58) 



Notice that for the integer flux vacuum we have A(oo) = 
1, but for the fractional flux vacuum ^4(oo) becomes frac- 
tional. 




FIG. 1: The D-type straight vortex with n = 1 and k = 1 
which has 4n/g flux. Three solutions are shown: a — /3 = 
7 = (solid lines), a = /3 = 0, 7 = 0.005 (dashed lines), and 
a = 7 = 0, /3 = —0.005 (dotted lines). Here the unit of the 
scale is 1/p and we have put X/g 2 = 2. 



At this point one might worry about the apparent sin- 
gularity at the core in the gauge potential when ^4(0) 7^ 0. 
But this singularity is a coordinate singularity which can 
easily be removed by a gauge transformation. Indeed one 
can always choose a gauge where ^4(0) becomes zero to 
remove the coordinate singularity. But notice that the 
gauge transformation also changes A(oo) by the same 
amount, leaving A(oo) — A(0) invariant. 

The existence of two types of boundary conditions in 
two-gap superconductor has an important impact. To 
understand this notice that the magnetic flux of vortex 
is given by 



$ = j> A^dx" = (A(oo) - A(0)) 



(59) 



Now, it is clear that the magnetic flux becomes fractional 
when A(oo) is fractional, which happens when < >=/= 1 
(or equivalently < | <^>2 1 >^ 0). As importantly, when 
A{po) = 1 the magnetic flux becomes 27r(n + k)/g with 
A(0) = —k/n. This was impossible in ordinary super- 
conductor. Now we classify the magnetic vortex in terms 
of the flux. 



A. 47r/(?-flux vortex 



Let us choose the Dirichlet boundary condition at the 
core and the integer flux vacuum at the infinity d, Il3| . 
With k = 1 we require 



p(0) = 0, p(0)#0, p(oo) 



' 2(p + 7 ) 
(A + a + f3) 
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FIG. 2: The D-type straight vortex with n = 1 and k = 2 
which has 6n/g flux. Three solutions are shown: a = j3 = 
7 = (solid lines), a = j3 — 0, 7 = 0.005 (dashed lines), and 
a = 7 = 0, /3 = —0.005 (dotted lines). Here the unit of the 
scale is 1/p and we have put X/g 2 — 2. 
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FIG. 3: The N-type straight vortex with n = 1 and 2-7r/<? 
flux. Three solutions are shown: a — (3 = 7 = (solid 
lines), a = /3 = 0, 7 = 0.005 (dashed lines), and a = 7 = 0, 
/3 = —0.005 (dotted lines). Here the unit of the scale is 1/p 
and we have put \/g 2 = 2. 



/(0) 
A(Q) = 



71", /(oo) = 0, 

1 



A(oo) = 1. 



(60) 



With this we can integrate I152H to find the vortex solu- 
tions. The solutions with n = 1 with different parameters 
are shown in Fig. ^ We call this a D-type vortex, be- 
cause this comes from the Dirichlet boundary condition 
at the core. Both tpi and 4>2 start from zero at the core. 
However, notice that cf>i approaches the finite vacuum 
value but 4>2 approaches zero at the infinity. So <f>2 has 
a maximum concentration at a finite distance from the 
core. This is a generic feature of a D-type vortex. 

With l|t)U|) the magnetic flux is given by 



= 1 



dnAmd X 



1 \ 27m 2tt , 

- — (n + 1 , 

n> g 9 



(61) 



so that when n = 1 the vortex has An/g flux. Moreover, 
the solution has a non-Abelian topology. To see this no- 
tice that n defines a mapping 7r 2 (S' 2 ) = n from the com- 
pactifed rry-plane S 2 to the CP 1 space S 12 . Clearly this 
is non-Abelian. 



In general we may require 



p(o) = d Ao) = 

ag 




(62) 



and obtain a different D-type vortex whose magnetic flux 
is given by 



1 



(n + k). 



(63) 



27m 27T 

.9 5 

The 67r/(7-flux vortex with n = 1 and /c = 2 is shown 
in Fig. |21 This tells that there exist infinitely many D- 
type vortices which have the same topology 7T2(S' 2 ) = n. 
Again this is completely unexpected. 



B. 27r/g-flux vortex 



Now we choose the Neumann boundary condition at 
the core and the integer flux vacuum at the infinity d, 
0, 



p(0)^0, p(0) = 0, p(oo) 



/(0) = 
A(0) 



TV, 

- 0, 



/(oo) 
A(oo) 



/ 2(p + 7 ) 
/ (A + « + /?)' 

0, 
= 1. 



(64) 



and find the vortex solutions. The solutions with n = 1 
but with different parameters are shown in Fig. [3J We 
call this a N-type vortex, because this comes from the 
Neumann boundary condition at the core. In this case 
(pi behavior is the same as before. But notice that <f>2 has 
a maximum concentration at the core, and approaches 
zero at the infinity. This is a generic feature of a N-type 
vortex. 

The magnetic flux of the vortex is given by 



$ = 



Fg^pd x — 



dgAipd x 



9 




magnetic flux is given by 



FIG. 4: The D-type straight vortices with n = 1 and k = 1 
which have a fractional flux with a = 7 = 0, /3 = 1.0 (solid 
lines), a = 0, /3 — X, 7 = —0.2 (dashed lines), and a — 
—0.25, j3 = A, 7 = (dotted lines). Here the unit of the 
scale is 1/p and we have put X/g 2 = 2. 



2nn 
9 



(65) 



so that it has the same flux as the Abrikosov vortex. But 
notice that the topology of the CP 1 field h is still non- 
Abelian as before, ^{S 2 ) = n. The reason why there 
exist two types of vortices which have different magnetic 
fluxes but have the same topology is that the magnetic 
flux is determined by the boundary condition A(oo) — 
A(0), not by the topology. The topology assures only the 
quantization of the flux, and does not determine what is 
the unit flux quantum. 



Fractional flux vortex 



This is possible when we have the fractional flux vac- 
uum at infinity 



P(°°)=^lW^, cos /( oo). 2 ^-^ 



4/3A - a 2 



2/3fi — off 



AM = lgM±M^to), (66) 

x ' 2 2/3^-aj y ' 

At the core we can impose either the Dirichlet condition 
(|55|l or the Neumann condition (|57l) . 

We consider two special cases: 
1. An = A 2 2 {a = 0). In this case we have 

p(oo) = y5 cos/(oo) = ^. (67) 
So with the Dirichlet boundary condition at the core the 



/ 7A 1 k \ 27m 
$ = — 1 1 — 

\2/3fi 2 nJ g ' 



(68) 



but with the Neumann boundary condition at the core 
the magnetic flux is given by 



\2Pfi 2) g 



(69) 



Clearly they are fractional. 

2. A12 = (A = 0). In this case two condensates (/)% and 
cj>2 have no direct coupling, and we have 



< \h\ >= Jg- (70) 




With this we have 



p(oo) 



l(2A-a)( M + 7) 



A(oo 



2 2^X - «7 



(71) 



So with the Dirichlet boundary condition at the core the 
magnetic flux is given by 



$ = 



■1 (2A-a)(M + 7 ) t k\2 
.2 2A/i — aj 



11/ g 



(72) 



but with the Neumann boundary condition at the core 
the magnetic flux is given by 



^_ l (2A-a)fc + 7) 2nn 
2 2\[i 



«7 



9 



(73) 



Again they are fractional, in spite of the fact that </>i and 
4>2 have no direct coupling. This is because they are cou- 
pled through the electromagnetic potential, which tells 
that the two-gap superconductor is not a naive superpo- 
sition of two one-gap superconductor. 

According to the different boundary condition at the 
core there are two types of fractional flux vortices, D-type 
and N-type. These fractional flux vortices with n — 1 are 
plotted in Fig. 0] and Fig.0 The fractional vortex is also 
topological, but the topology of the fractional vortex is 
different from that of integer flux vortex. Notice that 
for the fractional flux vortices the TT2(S 2 ) topology of n 
becomes trivial, ^(S* 2 ) = 0. This is because h does 
not cover the target space S 2 fully. But in this case 
we still have a £7(1) topology 7Ti(S' 1 ), the topology of 
the J7(l) symmetry which leaves h invariant. And this 
Abelian topology describes the topology of the fractional 
flux vortex. So the topology of the fractional flux vortices 
is the same as that of the Abrikosov vortex. 

An important feature of the fractional flux vortex is 
that the energy per unit length of the vortex is logarith- 
mically divergent, which can be shown from the Hamil- 
tonian l|51|). This is because the fractional flux vortex 
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FIG. 5: The N-type straight vortices with n = 1 which have 
a fractional flux with a — 7 = 0, f3 = 1.0 (solid lines), a = 
0, f3 = X, 7 = -0.2 (dashed lines), and a = -0.25, /3 = 
A, 7 = (dotted lines). Here the unit of the scale is 1/p and 
we have put X/g 2 — 2. 



does not mean that two-gap superconductor can not ad- 
mit an Abrikosov vortex. With / = tt (or / = 0) and 
a = (3 = 7 = 0, H52f> describes an Abrikosov vortex. This 
is because with (f>i — (or with (j) 2 = 0) the two-gap 
superconductor reduces to an ordinary superconductor. 



V. HELICAL VORTEX 

In this section we show that the above non- Abrikosov 
vortices can be twisted to form a twisted magnetic vortex. 
With the twisting we obtain the helical vortex which is 
periodic in z-coordinate. To show this we choose the 
following ansatz pl.lT3|. 



p{q), i = 



( f(o) , ■ \ 
cos „ exp(— imp) 



sinMexpR.) 



-4, 



-(nA 1 {g)d t ,ip + mkA 2 (g)d t ,zy (74) 



has a non-vanishing neutral current k^ at the infinity. 
This, however, does not make the fractional flux vortex 
unphysical. In laboratory setting one can observe such 
vortex because one has a natural cutoff parameter A fixed 
by the size of the superconductor, which can effectively 
make the energy of the fractional flux vortex finite. In- 
deed in He 3 superfluid one often encounters the vorticit; 
vortex whose energy is logarithmically divergent 



rorticity 

Emu. 



The existence of fractional flux vortices in two-gap 
superconductor has been pointed out before in London 
limit 15]. Our analysis in this paper shows that the 
London limit does not fully describe the vortex in two- 
gap superconductor. This is because the magnetic flux is 
determined by the boundary condition at the origin (as 
well as the boundary condition at the infinity). Clearly 
the existence of two types of vortices which have differ- 
ent core structure and different magnetic flux can not be 
understood in London limit. 

In this section we have shown that the two-gap su- 
perconductor can have totally different magnetic vor- 
tices which can not be found in ordinary superconductor. 
There are two types of vortices, D-type and N-type, and 
both have two different topologies, n 2 (S 2 ) and tt 1 {S 1 ), 
which describe the vortices. The integral flux vortex 
is described by the ir 2 (S 2 ) topology, but the fractional 
vortex is described by the 7Ti(S' 1 ) topology. As impor- 
tantly, there are infinitely many different vortices within 
the same topological sector. Moreover, the magnetic flux 
of the D-type vortex is larger than that of the N-type 
vortex by a factor 2irk/g, so that in the same topolog- 
ical sector the D-type vortex has more energy than the 
N-type vortex. 

Obviously all these vortices are non- Abrikosov. This 



Obviously the ansatz is periodic in z-coordinate, with the 
period 2t: jk. 

With the ansatz we have 

(sin /(g) cos(rnp + mkz) 
sin f(g) sm.(n<p + mkz) 
cos f{g) 

_ cos /(g) + l d m/; cos /(g) Zl 9 z 



= gp 



(n(A 



cos / + 1 . 

1 n 



ik(A 2 



cos / — 1 « 



2 ( (A f COS / + l \ Q 

\n\A-i_ cos / ) d M tp 

ik(A 2 cos / + — -)d^z 



gp 



and the following Hamiltonian 

^ = ^ 2 + ^ 2 (/ 2 + (^+- 2 fc 2 )sin 2 / 



(75) 



A x 



cos / + 1 \ 2 



, 2 fc 2 (A 2 -^ 1VJ 



+ ^k 2 A 2 



2^2 



+ ^ (p2 _^y cos/ + ^4 cos 2 / 



{p z -f) 

2 



p_ 

2A' 



(76) 



With this (|4*T|> becomes 



P+-P- 

e 



i(/ 2 + (^ + m 2 fc 2 )sin 2 /) 



cos / + 1 . 



-(A, - —L-^f + m 2 k 2 (A 2 - ~ y 



cos / — 1. 2 
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FIG. 6: The D-type helical vortex with 47r/<7-flux along the 
vortex with m = n = 1 (with k = 1). Three solutions are 
shown: a = j3 = 7 = (solid lines), a = /3 = and 7 = 0.005 
(dashed lines), a = 7 = and /3 = —0.005 (dotted lines). 
Here the unit of the scale is 1/p and we have put k — 0.12p 
and \/g 2 = 2. 



(P 2 



Pi) 



a ( 2 2 7 \ n , P 2 2 r 

A a A 



/+( 1 + 2^)/-2[^(A 1 -i) 
p p L g A 2 



( 



- + 5) 



sin/ 



1 



fc 2 (A 

27 ~(f +/5cos/)p 2 )sin/ 
cos / + 1 



<?V(a 



0. 
0. 



(77) 



A a - -At 
Q 

Y 1 : 99/, cos f — 1 

a 2 + - a 2 - <?y a 2 — j - — 

g \ 2 

This is an obvious generalization of (|52[1 . 

To obtain the helical vortex we first consider the in- 
teger flux boundary condition at the infinity 



p(oo) 



2(M + 7) 



y (A +0 + a)' 
Ai(oo) = 1, A 2 (oo) 



/(«>) 
0. 



(78) 



Just like the straight vortex, there are two types of 
boundary conditions at the core. Here we consider only 
the case n = 1 for simplicity: 
A. Dirichlet boundary condition 



P (0) = 0, /(0) = 7T, 

Ai(0) = -1, A 2 (0) = 0. 
B. Neumann boundary condition 

p(0) = 0, /(0) = 7T, 

Ai(0) = 0, A 2 (0) = 0. 



(79) 



(80) 



FIG. 7: The N-type helical vortex with 27r/g-flux along the 
vortex with m = n = 1. Three solutions are shown: a = (3 = 
7 = (solid lines), a — (3 = and 7 = 0.005 (dashed lines), 
a = 7 = and /3 = —0.005 (dotted lines). Here the unit of 
the scale is 1/p and we have put k = 0.12p and X/g 2 = 2. 



With the Dirichlet boundary condition we have the D- 
type helical vortex which has 477/17- flux along the vor- 
tex shown in Fig. but with the Neumann boundary 
condition we have the N-type helical vortex which has 
27r/<7-nux along the vortex shown in Fig. \7\ 

For the fractional flux helical vortex we impose the 
fractional flux boundary condition at the infinity 



p(oo) = 2 



2/3fi — off 
4/3 A - a 2 ' 



cos/(oo) = 



27A — afi 

2(3 jj, — ' 



1 2( 7 A + /3^)-q(m + 7) 
Ai 00 = , 

V ; 2 20n-ary 
, . , 1 2( 7 A - (3p) - qfc - 7) 

Now, with the Dirichlet boundary condition at the core, 
we obtain the D-type helical vortex which has a frac- 
tional flux along the vortex shown in Fig. |SJ But with 
the Neumann boundary condition (|80|l at the core, we 
obtain the N-type helical vortex which has a fractional 
flux along the vortex shown in Fig. EI Notice that, just as 
the fractional flux straight vortex, the energy per one pe- 
riod of the fractional flux helical vortex is logarithmically 
divergent. 

A new feature of the helical vortex is that the mag- 
netic flux becomes helical. Indeed the ansatz (|74(l tells 
that the magnetic flux can be decomposed to the one 
along the vortex and the other around the vortex 0, 0] 



-An 



J? - nAl TP -H^ 

r gip ? -t'zg 

9 Q 9 
so that we have two magnetic fluxes linked together, 

27rn 
9 



(81) 



**= J F^gdgdip = (Ai(oo) - Ai(0) 
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FIG. 8: The D-type helical vortices which have a fractional 
flux with m — n = 1 (with k = 1). Two solutions with 
a — 7 = 0, f3 — 0.5 (solid lines), and a = 7 = 0, f3 = 1.0 
(dashed lines) are shown. Here the unit of the scale is 1/p 
and we have put \/g 2 = 2, k — O.lp. 



FIG. 9: The N-type helical vortices which have a fractional 
flux with m = n — 1. Two solutions with a = 7 = 0, /? = 0.5 
(solid lines), and a = 7 = 0, /3 = 1.0 (dashed lines) are 
shown. Here the unit of the scale is 1/p and we have put 
X/g 2 = 2, k = O.lp. 



J F se dzdQ = -(A 2 (o O )-A 2 (0))^. (82) 

Obviously $^ is due to the helical structure of the vortex, 
which becomes fractional in general. 

Another important feature of the helical vortex is 
that the electromagnetic current which is responsible 
for the Meissner effect also becomes helical. In particu- 
lar, it has a non-trivial electromagnetic current jg along 
the vortex which generates the magnetic flux <I>^, in ad- 
dition to the usual electromagnetic current jq> around 
the vortex which is responsible for $£. But notice that 
the total electromagnetic current i% along the vortex be- 
comes zero. This, together with (|50f) . tells that <t>\ and 

cj)2 generate non-vanishing electromagnetic currents 

(2) 

and v~ which flow oppositely and cancel each other. In 
this sense we may call the helical vortex superconduct- 
ing, even though it has no net electromagnetic current ij 
along the vortex 0, . 



VI. MAGNETIC KNOT IN TWO-GAP 
SUPERCONDUCTOR 

Clearly the helical vortex is unstable unless the pe- 
riodicity condition is enforced by hand. Nevertheless it 
has an important implication, because the helical vortex 
predicts the existence of a topological knot in two-gap su- 
perconductor. This is because we can make it a twisted 
magnetic vortex ring smoothly bending and connecting 
two periodic ends together. The resulting twisted mag- 
netic vortex ring becomes a knot whose topology is de- 
scribed by the Chern-Simon index of the electromagnetic 



potential [113. 

There have been two objections against the existence 
of a stable magnetic vortex ring in Abelian superconduc- 
tor. First, it is supposed to be unstable due to the tension 
created by the ring [22j . Indeed if one constructs a vor- 
tex ring from an Abrikosov vortex, it becomes unstable 
because of the tension. But we can easily overcome this 
difficulty by twisting the magnetic vortex first and con- 
necting the periodic ends together. In this case the non- 
trivial twist of the magnetic field forbids the untwisting of 
the vortex ring by any smooth deformation of field config- 
uration, and the vortex ring becomes a stable knot. The 
other objection is that the Abelian gauge theory is sup- 
posed to have no non-trivial knot topology which allows 
a stable vortex ring. This again is a common misconcep- 
tion. As we have seen, the theory has a well-defined knot 
topology TT3(S 2 ) described by the Chern-Simon index of 
the electromagnetic potential. This tells that there is no 
reason whatsoever why the Abelian superconductor can 
not have a topological knot. 

To demonstrate the existence of a topological knot in 
the two-gap superconductor, we introduce the toroidal 
coordinates (rj,j,(p) defined by 

x = — sinh r\ cos (f, y = — sinh i] sin ip, 
a 

Z = — sin 7, 

D = cosh?; — cos 7, 
2 

da 2 = ^ (drj 2 + d-y 2 + sinh 2 rjdip 2 ^ , 
a 3 

d 3 x — — 5- sinh rjdr]d'yd(p, (83) 
where a is the radius of the knot defined by 77 = 00. 
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Notice that in toroidal coordinates, rj — 7 = represents 
spatial infinity of R 3 , and 77 = 00 describes the torus 
center. 

Now we choose the following ansatz, 
/ (77 j 7) 

I / cos ■ — exp(— im</?) \ 

v \sin — - — exp(zno;(?7, 7)) / 

^ = -A (?7, 77 + -Ai (77,7)^7 



+—M{r],l)d^ip. 
9 



With this we have 

(sin f(rj, 7) cos(nw + imp) \ 
sin/(?7, 7) sin(ncj + rayj) , 
COS/(T7,7) / 
777 77 



(84) 



777- 777 

F lv = —d 1 A 2 , F lpr , = - — d n A 2 . (85) 
Notice that, in the orthonormal frame (77,7, 0), we have 



D 



D 



D 



A, — _A A. — A A „ — 4 

a a a sinh 77 

77Z? 2 

— — , a 7 A 2 , 
gar sinn 77 

mfl 2 

For? — 9 ■ ; OnAz- 

ga smn 77 

Next, we adopt the SU(2) symmetric potential with 
An = A22 = A12 = A and //1 = /j,2 = // for simplicity. 
In this case we have the following knot equation 



^ + ^ + Vsinh7 ? D £> ^ 



+ sin 2 /(n 2 (9^) 2 + n 2 (a^) 2 + -^-) 



n 2 Uo 



m cos / — 1 n '-' 

sinh 2 77 ^ 2 



cos / + 1 
2 

A a 2 



(a t) /) 2 + (9 7 /) 2 
a A 2 , 2/, cos/ + l 

c^w) +77 lAi 



2 ^ 



2 sin / n 2 (Aq - ^d v u?j 8 v lu + n 2 (a x - ^d y u?j d y uj + ^™ 2 ( A ? + 



rl 2 , rl 2 , /£oshj7 _ smh7? \ _ sin7 
^ + ^ + Vsinh77 D r" D 7 



+^ 7 ($,/4> + a Y M l ) =0, 
sin; V / 



a 



S1117 



cosh 77 sinh 77 
sinh 7/ 



— y ^a 7 a q - a^Ai j = 752 5 2 p 2 ( A) - 

2 



cos/ + l 
2 9 >7 w j> 



D 



cos / + 1 



a 7 w , 



2 2 /cosh?7 sinh77\ sin7 



'7J ' "7 



4 = ^2.9 2 P 2 (A 2 



cos / - 1 



(86) 



Moreover, from (JIJ and l|84|l we have the following Hamiltonian for the knot 



n 



LP 
2a 2 



(d v p) 2 + (d,p) 2 + ip 2 ((5 r) /) 2 + (d 7 /) 2 ) + ip 2 sin 2 /(77 2 (a^) 2 + n\d^f 



sinh" 77 
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2a 2 



2 COS/ + 1 2 2 C0S/ + 1 2 TO 2 COS / - 1 2 

n [A d r ,uj) +n{Ai ; • ' •■ 1 1 



COS/+ 1 x 2 , 171 , A 

4 smn ry 



D 4 



2. 9 2 a 4 



2 

n 2 (a,Ai - 9 7 A ) 2 + -^2- (W2) 2 + (9 7 A 2 ) 2 ) 
sinn v 



+ ^(P 2 -P 2 ) 2 - 



(87) 



Minimizing the energy we reproduce the knot equation 

©. 

In toroidal coordinates, 77 = 7 = represents spatial 
infinity of R 3 , and 77 = 00 describes the torus center. So 
we can impose the following Neumann boundary condi- 
tion 



p(0,0) = p, 
/(0,7) = 0, 
w(*7,0)=0, 
A (0, 7 ) = 0, 
Ai(0 >7 ) = 1, 
A 3 (0, 7 ) = 0, 



p(oo,7) = 0, 

/(00,7) = 7T, 
w(?7, 27r) = 27T, 

4,(00,7) = 0, 
Ai(oo,7) = 0, 
4(00,7) = -1, 



(88) 



to obtain the desired knot. A numerical integration of 
(|86|l with the boundary conditions l|88|) is difficult to per- 
form. But we can obtain the actual knot profile of p, f, 
and u) by minimizing the energy. From (|87|) the knot 
energy is given by 




FIG. 10: Fhe energy profile of a N-type knot with m = n = 1. 
Here we have put \/g 2 = 2. 



E = J TLjj^ smhrjdrjd'ydip. 



(89) 



We find that, for m = n = 1, the radius of knot which 
minimizes the energy (|89|) is given by 



1.2 



2/i 



(90) 



From this we obtain the three-dimensional energy profile 
of the lightest axially symmetric knot in two-gap super- 
conductor, which is shown in Fig. 1101 

With this we can estimate the energy of the axially 
symmetric knot, 



The flux is quantized in the unit of 2ir/g, but this is 
due to the SU (2) symmetric potential and the Neumann 
boundary condition (|88|l . In general they can be frac- 
tional. But independent of this the two fluxes are linked, 
whose linking number becomes mn. This is important. 

This confirms that the knot can be viewed as a twisted 
magnetic vortex ring, where the linking of two magnetic 
fluxes provides the knot topology. There is a natural 
candidate which can describe this topology of twisted 
magnetic vortex ring, the Chern-Simon index of the elec- 
tromagnetic potential which describes the tts(S 2 ) topol- 
ogy. We can calculate the knot quantum number from 
the Chern-Simon index, and find 



E ~ 51 



(91) 



We can also calculate the magnetic flux of the knot. Since 
the flux is helical, we have two fluxes, the flux <I> 7 passing 
through the knot disk of radius a in the ccy-plane and the 
flux <I>^ which surrounds it. From the knot solution we 
find 



$ 7 = 



F~ fv d'ydtp = 0. 



F vv dipdr] = 



F^dndj 



2m7r 
9 

2mr 



(92) 



Qcs 



mn 
8^2 



s^AiFjkdx 3 = 



(93) 



This confirms that the Chern-Simon index is indeed given 
by the linking number of two magnetic fluxes $ 7 and $^ . 

It has been asserted that the knot topology is de- 
scribed by the tt 3 (S 2 ) topology of the CP 1 field n 0- 
We emphasize that this is only partially true, which be- 
comes correct only when the knot carries an integer mag- 
netic flux. Indeed, in this case h acquires a non-trivial 
knot topology n 3 (S 2 ), which is given by @ 



mn 



(drjfdyU) — d-yfd v uj) sin fdrjd'ydcp 
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/ sax t at du = mn, 

4tt 



(94) 



where the last equality comes from the boundary condi- 
tion Notice, however, that the Kz{S 2 ) topology of h 
becomes trivial when the knot has a fractional flux. This 
is because the fractional flux knot comes from a frac- 
tional flux vortex, which has a trivial ^(S 12 ) topology. 
This means that the ^(S 2 ) topology of n can not de- 
scribe the topology of a fractional flux knot. In contrast, 
the Chern-Simon index of the electromagnetic potential 
can still provide the non-trivial tt^{S 2 ), even when the 
knot has a fractional flux. This is because the Chern- 
Simon index describes the knot topology of the twisted 
magnetic flux. 



VII. JOSEPHSON INTERACTION 



It has been well-known that two-gap superconductor 
may allow the prototype Josephson interaction |14| 



77(0102 + 0102); 



(95) 



where rj is a coupling constant. But we can consider a 
more general quartic Josephson interaction, 



Vj = fj (<pl(j)2 exp(i6»i) + $,01 exp(-i0i)) 
+fh [(01 fa? expW + (020i) 2 exp(-ifla) 
-^ 2 (|^i | 2 + |0 2 | 2 ) [{4>\(P2 exp(i0 3 ) + 05&esp(-i0 3 )) 
+%(l0i| 2 -|02| 2 )[(0i02cxp^ 4 ) 



+0201 exp(— 16*4)) 



(96) 



Clearly the Josephson interaction breaks the U(l)xU(l) 
symmetry of the potential JSJ down to U(l). In general 
it is not easy to accommodate this type of generalized 
Josephson interaction. But for simpler Josephson inter- 
actions we can accommodate them within the framework 
of the potential JSJ). To understand how, notice that the 
potential (J5J already has a Josephson interaction in the 
sense that it allows an interband transition between 
and 4>2 when A12 is not zero. This implies that the above 
Josephson interaction could be included in the A12 inter- 
action. 

To exploit this point we generalize the U(l) x U(l) 
symmetric potential (JSJ to 



v -i-y = v + v u 

Vi=V (0102 exp(«0) + 020i exp(-id)) 



+?7i 



(05; 2 ) 2 exp(2i0) + (0^0i) 2 exp(-2*0)) 



+?72(|0i| 2 + |0 2 | 2 ) (0i0 2 exp(^) + 0^01 exp(-»0)) 



+%(|0i| 2 -|0 2 | 2 ) (0102 exp^) 



+020i exp(-i0)) 



(97) 



and introduce a new doublet ip with an SU (2) transfor- 
mation of 6, 



a 



a 



M = 



cos — exp(z6) — sin — exp(ic) 
sin — exp(— ic) cos — exp(— ib) 



tana 



V 

1 



Now, we can show that when 
9 = -b 



<h 



2 7 2 _ v 2 



2^2 — jy 2 
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(99) 



the potential l|97[) can be written as 



Aii, 
2 



^^ir+A'^iiw 



-MilV'il 2 -^1^2 



A' 



{\M 2 + \H 



y ) 2 + f(l^, 



+§(lv>ll 2 



l^2| 2 ) 2 - 

2A' 



(100) 



where 



A' = A 



/i = fX, 



0' = 



2(7 2 



2^ 2 — fj 



7' = \ll 2 + V 2 - 



(101) 



So in terms of ipi and -02 the potential 19 7|) becomes 
the potential (|35|l which has no Josephson interaction. 
This means that, with (|99J) . we can formally absorb the 
Josephson interaction to A12 interaction. From this we 
concludes that the presence of the Josephson interaction 
does not affect the existence of the topological objects in 
two-gap superconductor. 

This does not mean that the Josephson interaction 
does not affect the topological solutions. On the contrary, 
it does change the shape of the solutions drastically. This 
is because under the transformation l|98|l the profile of 
0i and 02 change drastically. To demonstrate this we let 
a — (3 = for simplicity, and adopt the potential which 
has the following Josephson interaction 



V 



^(I0l| 2 + \<P2\ - -I )\\"\\ 

+ry(0J0 2 exp(i0) + ^ exp(-i0)). (102) 



102 



V\2 

X 1 



7(l0l| 2 -|02| 2 ) 
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(a) \<t>i\ 2 



(b) |02| 2 



FIG. 11: The density profile of |</>i| 2 and \4>2^ of the N-type magnetic vortex in the presence of Josephson interaction. Here 
we have put p = 1, 7 = 0.05, 77 = 0.25, and \/g 2 = 2. 



Notice that, in terms of ip, the potential is written as 



2 _ M\2 

Y 



(103) 



so that it clearly has two types of straight vortex solution 
of the following form 



„ / cos — exp(— imp) " 



. / 



p = p(q), 



n 



(104) 



= -A(g)d^. 
Now, in terms of <fi, the solution acquires the form 



c fa,. .,s 

4i = cos — cos — exp(— imp — ib) 

f a 
+ sin — sin — exp(i9 + ib), 

^2 = — cos — sin — exp(— imp — id — ib) 

f Q 

+ sin — cos — exp(i&). (105) 

So, in the N-type vortex both 4>i and <^2 have a non- 
vanishing concentration of at the core in the presence 
of the Josephson interaction. The density profile of cf>± 



and (f>2 of the N-type 27r/g-flux vortex with 6 = and 
6 = is plotted in Fig. ^2 Obviously the solution is not 
axially symmetric. More importantly the vortex appears 
as a "bound state" of two vortices made of <pi and 4>2- 
This confirms that the Josephson interaction does not 
prevent the the existence of two types of magnetic vortex, 
but changes the profile of the solutions drastically. We 
notice that a similar vortex has been discussed in two- 
component Bose-Einstein condensate |23j. 

Furthermore we can construct a helical vortex by 
twisting the above vortex and making it periodic in z- 
coordinate. In this case the helical vortex becomes a 
"braided" magnetic vortex made of two vortices of <f>\ 
and <f>2- To see this we let 



( f r • \ 
cos — exp {—imp) 

f 

. sin — exp {imkz) 



= -A 1 {g)d^+ —A 2 {g)d f _ l z, 
9 9 



(106) 



and again obtain two types of vortex. In this case the 
solution has the following particle densities for <pi and 



i*i 2 



4 V 
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:COS/ 



\Jl 2 + V 



^ sin / cos(niy9 + mkz + 9 + 2b)^j , 
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71 sin / cosfny + mfci + fl + 26)^) (107) 

vV + ?? 2 y 

Clearly this shows that in the presence of the Josephson 
interaction the helical magnetic vortex becomes a braided 
vortex in which 0i-flux and 02-flux are braided together. 

Now, it goes without saying that we can make a 
braided knot, a twisted vortex ring, with the braided 
magnetic vortex. This tells that the Josephson interac- 
tion makes the topological objects in two-gap supercon- 
ductor more interesting. 



VIII. NON-ABELIAN SUPERCONDUCTOR 

So far we have discussed an Abelian gauge theory of 
two-gap superconductor. But our analysis implies that 
the doublet (</>i, 02 ) can be treated as an SU(2) doublet. 
Indeed, when a — (3 — 7 — 0, the Lagrangian (|37|l has 
an exact SU(2) symmetry. Even when there is no SU{2) 
symmetry, one may still regard that the theory has an 
approximate SU (2) symmetry which is broken by the a, 
/3, and 7 terms. In this sense one may conclude that the 
Abelian two-gap superconductor has a (broken) SU{2) 
symmetry. On the other hand the Ginzburg-Landau La- 
grangian H37fl is still based on the Abelian electromag- 
netic interaction. This leads us to wonder whether one 
can have a genuine non- Abelian superconductor in which 
the superconductivity is described by a non- Abelian dy- 
namics. 

To discuss this issue, notice that in the above two- 
gap superconductor the two condensates <pi and 02 carry 
the same charge, because the doublet is coupled to the 
Abelian electromagnetic field. Now we show that when 
the two condensates are made of opposite charges (made 
of one electron-electron pair condensate and one hole-hole 
pair condensate) the two-gap superconductor can be de- 
scribed by a genuine non- Abelian SU (2) gauge theory. 
Moreover, we show that this type of non-Abelian super- 
conductor also allows a non-Abrikosov magnetic vortex 
and topological knot identical to what we have discussed 
in this paper. 

To construct a theory of non-Abelian superconduc- 
tivity which is based on a genuine non-Abelian gauge 
theory, we need to understand the mathematical struc- 
ture of the non-Abelian gauge theory. In non-Abelian 
gauge theory one can always decompose the gauge po- 
tential into the restricted potential and the valence 
potential A M . Consider the SU(2) gauge theory and let h 
be a gauge covariant unit triplet which selects the charge 
direction of SU(2). In this case we have the following 



decomposition [TsUI^. 

A,j, = A^n n x dfj,n + X^ = A^ + A M , 

(A„ =h-A^ h 2 = 1, h-X„ = 0), (108) 

where A^ is the "electric" potential. Notice that the 
restricted potential is precisely the potential which leaves 
h invariant under the parallel transport, 

D^n = df^h + gA^ x h = 0. (109) 

Under the infinitesimal gauge transformation 

Sh = —a x n , SA^ = -D^a, (HO) 

one has 

5Au = -h ■ dud, SAn = -D^d, 
9 9 

SXft = -a X AV (111) 

This tells three things. First, A^ by itself describes an 
SU(2) connection which enjoys the full SU(2) gauge de- 
grees of freedom. Secondly, the valence potential A M 
forms a gauge covariant vector field under the gauge 
transformation. Most importantly, this tells that the 
decomposition is gauge-independent. Once the gauge co- 
variant topological field n is given, the decomposition fol- 
lows automatically independent of the choice of a gauge 

MM 

The importance of the decomposition (|108fl for our 
purpose is that one can construct a non-Abelian gauge 
theory, a restricted gauge theory which has a full non- 
Abelian gauge degrees of freedom, with the restricted po- 
tential A^ alone [3 ED ■ This is because the valence po- 
tential can be treated as a gauge covariant source, so 
that one can exclude it from the theory without compro- 
mising the gauge invariance. Indeed it is this restricted 
gauge theory which describes the non-Abelian gauge the- 
ory of superconductivity 0. 

Remarkably the restricted potential A^ retains all the 
essential topological characteristics of the original non- 
Abelian potential. In fact, h defines the ^(.S 2 ) topol- 
ogy which describes the non-Abelian monopoles and the 
tts(S 3 ) topology which characterizes the topologically 
distinct vacua 0, 0, 0|. Furthermore it has a dual 
structure, 

F^u = d^Av — 0,/Afi, 
H^v = --ft ■ (df.n x d v h) = d l ,C v - d v C^, (112) 

where C M is the "magnetic" potential. Notice that this 
is exactly the potential C M that we have introduced in 
(|40|l . This is an indication that the Ginzburg-Landau 
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theory of two-gap superconductor is closely related to 
the restricted SU (2) gauge theory. 

With these preliminaries we now demonstrate a non- 
Abelian superconductivity and non-Abelian Meissner ef- 
fect. Consider a SU(2) gauge theory described by the 
Lagrangian in which a doublet $ couples to the restricted 
SU{2) gauge potential, 



1 



C = -\b ll *\*-V{§,&)- 7 F. 



(113) 



The equation of motion of the Lagrangian is given by 

dV 



D^F^v = g 



L> 2 $ = — 

(D v $)iZ:& (114) 
2,1 2i J 



Let £ and r\ be two doublets which form an orthonormal 
basis, 



i\i = 1, rprj = 1, - tft = 0, 



ry or\ = —n, 



(115) 



(n-a) £ = (h-a) rj = -77, 
and let 

$ = (/)+e + 0-ry, (</»+ = 0_ = (116) 

With this we have the identity 



— — (C^h + -h x d^h) ■ a 



1 



2i 



&u + 7r(Cu« n x 9„n) • a 

2i g 



r, = 0, 



and find 



where 



(117) 



(118) 



A, 



2i 



9 



9 



From this we can express (|113ll as 
£ — _\n a . \2 _ \ n a, 12 



where 



-|^+| 2 -|^_| 2 -y(0 + ,0_) 
^(^0 + + t _0_) 2 -^ 2 ,, 



T = B A — 8 A 



(119) 



This tells that the restricted SU(2) gauge theory (|113|l is 
reduced to an Abelian gauge theory coupled to oppositely 



charged scalar fields d) + and <p- . We emphasize that this 
Abelianization is achieved without any gauge fixing. 

The Abelianization assures that the non-Abelian 
Ginzburg-Landau theory is not different from the Abelian 
Ginzburg-Landau theory. Indeed with 



(120) 



i^ 2 

4 



we can express the Lagrangian (|119|) as 
£ = -\D^ X \ 2 +M 2 X t X-^(x t *) 2 



(121) 



This is formally identical to the Lagrangian (|37[) of the 
Abelian two-gap superconductor discussed in Section III. 
The only difference is that here 4> and are replaced 
by x and A.u. This establishes that, with the proper 
redefinition of field variables l|lltj|) and l|118|) . our non- 
Abelian restricted gauge theory (|113fl can in fact be made 
identical to the Abelian gauge theory of two-gap super- 
conductor. This proves the existence of non-Abelian su- 
perconductors made of the doublet consisting of oppo- 
sitely charged condensates 0. As importantly our anal- 
ysis tells that the two-gap Abelian superconductor has 
a hidden non-Abelian gauge symmetry, because it can 
be transformed to the non-Abelian restricted gauge the- 
ory. This implies that the underlying dynamics of the 
two-gap superconductor is indeed the non-Abelian gauge 
symmetry. In the non-Abelian formalism it is explicit. 
But in the Abelian formalism it is hidden, where the full 
non-Abelian gauge symmetry only becomes transparent 
when one embeds the doublet (4>i,4>2) properly into the 
non-Abelian symmetry. 

Once the equivalence of two Lagrangians (|37|1 and 
(11131 is established, it must be evident that the non- 
Abelian gauge theory of two-gap superconductor also ad- 
mits a non-Abrikosov vortex and magnetic knot. This 
confirms the existence of a non-Abelian Meissner effect 
and non-Abelian superconductivity. All the above results 
of Abelian superconductor become equally valid here. 



IX. DISCUSSION 



In this paper we have shown that the two-gap super- 
conductor can admit non-Abrikosov vortex and topolog- 
ical knot. There are two types of non-Abrikosov vortex, 
D-type and N-type. The D-type has no concentration of 
the condensate at the core, but the N-type has a non- 
trivial profile of the condensate at the core. In terms of 
topology there are two, the non-Abelian topology ir^S 2 ) 
defined by h and the Abelian topology ni(S 1 ) defined 
by the invariant subgroup of h. And both D-type and 
N-type vortices exist within the same topological sector. 
In particular, we have infinitely many D-type vortices 
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which have the same topology. The magnetic flux of the 
vortices can be integral or fractional. The N-type vor- 
tex can have a 27rn/g-flux or a fraction of this flux, but 
the D-type vortex has 2?tk/g more flux than the N-type. 
And we have shown that these non-Abrikosov vortices 
can be twisted to form a helical vortex which is periodic 
in z-coordinate. 

Perhaps a most interesting topological object in two- 
gap superconductor is the magnetic knot, a twisted mag- 
netic vortex ring made of helical vortex. Our analysis 
suggests that we have two types of knot, the D-type and 
the N-type. They are made of two magnetic fluxes linked 
together, one flux along the knot axis and one flux along 
the knot tube. And the linking number of two fluxes 
provides the knot topology ^(S 12 ), which is described by 
the Chern-Simon index of the electromagnetic potential. 
The knot is stable dynamically as well as topologically. 
The topological stability follows from the fact that two 
flux rings linked together can not be separated by any 
continuous deformation of the field configuration. The 
dynamical stability follows from the fact that the flux 
trapped inside of the knot ring can not be squeezed out, 
which means that it provides a repulsive force against the 
collapse of the knot. Another way to understand this dy- 
namical stability is to notice that the supercurrent along 
the knot generates a net angular momentum around the 
knot axis. And this provides the centrifugal repulsive 
force preventing the knot to collapse. This makes the 
knot dynamically stable. 

The Josephson interaction makes these topological 
objects more interesting. The straight vortex becomes 
a bound state of two magnetic vortices made of two con- 
densates cf>i and 4>2, and the helical vortex becomes a 
braided magnetic vortex of two condensates. Moreover 
the knot acquires the form of a braided magnetic vortex 
ring. And we have two of them. 

It must be emphasized, however, that the actual mag- 
netic flux of vortex and knot is determined by the two- 
gap superconductor at hand because it is fixed by the 
parameters of the potential which characterizes the su- 
perconductor. Independent of this all two-gap supercon- 
ductors have two types of vortex and knot. On the other 
hand one must keep the followings in mind. First, com- 
pared with the N-type vortex the D-type vortex has more 
energy in general because the D-type carries more flux. 
This opens the possibility that, within the same topologi- 
cal sector, the D-type vortices could decay to the N-type 
vortices. Secondly, the energy (per unit length) of the 
fractional flux vortex and knot is logarithmically diver- 
gent, so that they can exist only when there is a cut- 
off parameter which makes the energy finite. This tells 
that the N-type 2n/g vortex forms the true finite energy 
ground state vortex of two-gap superconductor. 

Another important lesson from our analysis is that the 
non-Abelian dynamics could play a crucial role in con- 
densed matter physics. Indeed we have shown that we 



can actually construct an SU{2) gauge theory of super- 
conductivity which is mathematically equivalent to the 
Abelian gauge theory of two-gap superconductor. This 
means that, implicitly or explicitly, the underlying dy- 
namics of multi-gap superconductor can ultimately be 
related to a non-Abelian gauge theory. 

In this paper we have studied the topological objects 
in two-gap superconductor. But from our discussion it 
must be clear that similar topological objects should also 
exist in multi-gap superconductor in general. This is 
because the multi-gap superconductor is described by a 
multi-component condensate, which naturally accommo- 
dates the non-trivial non-Abelian topology. 

Clearly the above theory of two-gap superconduc- 
tor is closely related to the Gross-Pitaevskii theory of 
two-component Bose-Einstein condensate (BEC), which 
tells that similar topological objects can also exist in 
two-component BEC 0, ITTl Il2j . This is because in 
the absence of the electromagnetic interaction the above 
Ginzburg-Landau Lagrangian reduces to the Gross- 
Pitaevskii Lagrangian of two-component BEC. But there 
is an important difference. In two-component BEC only 
the N-type vortex and knot exist, because it allows only 
the N-type boundary condition |g, u3 • In this sense it 
is really remarkable that two-gap superconductor allows 
two types of topological objects. 

We close with the following remarks: 

1. Recently similar non-Abelian vortices and knots have 
been asserted to exist almost everywhere, in atomic 
physics in two-component BEC 0, 0, ^| , in condensed 
matter physics in multi-gap superconductors 0, UH > m 
nuclear physics in Skyrme theory 0, l2l| , in high energy 
physics in QCD [24|. The major difference here is that 
our vortex and knot are made of a real magnetic flux. We 
emphasize that at the center of these topological objects 
lies the baby skyrmion and the Faddeev-Niemi knot in 
Skyrme theory. In fact, one can show that our magnetic 
vortex and knot (as well as those in two-component BEC) 
are a straightforward generalization of the baby skyrmion 
and the Faddeev-Niemi knot 0, 0, . This is because 
both the Ginzburg-Landau theory of two-gap supercon- 
ductor and the Skyrme theory are described by a CP 1 
field n which obeys the same non-linear dynamics. 

2. From our analysis there should be no doubt that the 
non-Abrikosov vortex and the magnetic knot must exist 
in two-gap superconductor. If so, the challenge now is 
to verify the existence of these topological objects exper- 
imentally. Constructing the knot might not be a sim- 
ple task at present moment. But the construction of 
the non-Abrikosov vortex could be rather straightfor- 
ward (at least in principle) |25|. To identify the non- 
Abrikosov vortex, there are two points one has to keep in 
mind. First, the magnetic flux of the non-Abrikosov vor- 
tex need not be 2ir/g, and can be fractional in general. 
More importantly, there are two types of vortex, the D- 
type which has no concentration of the condensate at the 
core and the N-type which has a non-trivial concentra- 
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tion of the condensate at the core. These are the crucial 
points which distinguish the non-Abrikosov vortex from 
the Abrikosov vortex. With this in mind, one should be 
able to construct and identify the non-Abrikosov vortex 
in two-gap superconductor without much difficulty. 
3. The non-Abelian gauge theory of superconductivity is 
not just an academic curiosity. There is an excellent ex- 
ample of non-Abelian two- gap superconductor, the liquid 
metallic hydrogen (LMH) |2(|. Under high pressure the 
LMH becomes a superconducting state in low tempera- 
ture, due to the electron Cooper pairs. But in a lower 
temperature the proton Cooper pairs can coexist with 
the electron pairs. And obviously it has no Josephson in- 
teraction and probably a weak or no A12 interaction. So 
the LMH becomes an excellent candidate of non-Abelian 
two-gap superconductor. This implies that the LMH can 



have all the topological objects we have discussed in this 
paper. In particular it must have two types of vortex and 
knot. 

In this paper we have discussed the topological objects 
which we obtain with the ansatz (@^|, ifTljl. (|%4")l . or (|105jl . 

But we emphasize that there are other topological objects 
which can be obtained with different ansatz. These ob- 
jects and the physical implications of these topological 
objects in two-gap superconductor will be discussed in 
an accompanying paper |27| . 
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